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Abstract. We present a comprehensive X-ray study of four years of pointed .RJfT'iJ observations of GRS 1915+105 
in the x-state. We interpret the behavior of the hard power law tail spectrum as coming from inverse Compton 
scattering of soft disk photons on a thermally dominated hybrid corona above the accretion disk. GRS 1915+105 
shows a strong, variable reflection amplitude. As in other BHC and in Seyfert galaxies, a correlation between the 
power law slope and the reflection was found. Also, the radio fluxes at 2.25 GHz and 15 GHz correlate with the 
power law slope, thus revealing a connection between the outflowing matter and the comptonizing region in the 
X-states. 
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1. Introduction 

With the discovery of superluminal ejections (Mirabel & 
Rodriguez 1994) from the galactic transient source GRS 
1915+105 and GRO J1655-40, much observational and 
theoretical attention has been directed to the field of mi- 
croquasars. Microquasars are thought to be downscaled 
analogs to quasars exhibiting much smaller time scales 
and are therefore potential laboratories for studying ac- 
cretion and relativistic jets near black holes (Mirabel et al. 
1992). The most prominent object of this type is the galac- 
tic X-ray binary system GRS 1915+105 which shows dra- 
matic variability (Greiner et al. 1996) in light curve, quasi- 
periodic oscillations, phase lags and coherence behavior 
(Morgan etal. 1997, Muno etal. 2001). GRS 1915+105 
is the most energetic object known in our galaxy with a 
luminosity of ^10^^ erg/s in the low state. 

GRS 1915+105 harbors a black hole of 14 Mg (Greiner 
etal. 2001b) making it the most massive stellar black hole 
known. Being in the galactic plane and at a distance of 
~12 kpc the source suffers a very high extinction in the 
optical band of the order of 25-30 mag (Greiner et al. 1994, 
Chaty et al. 1996). Greiner et al. (2001a) found the donor 
to be an K-M III late type giant as determined from ab- 
sorption line measurements in the near infrared. 
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Black hole transients generally exhibit different states 
of intensity and spectrum: the high/soft state with a 
prominent disk component and a weak or negligible steep 
power law tail, the low/hard state with negligible accre- 
tion disk and flat power law. In the intermediate and very 
high state both the accretion disk and a power law are 
seen. 

Originally discovered by Granat (Castro- Tirado et al. 
1992), GRS 1915+105 was extensively monitored by the 
RXTE since 1996 and a number of investigations of these 
data have been published (e. g. Morgan et al. 1997, Belloni 
etal. 1997 & 2000, Muno etal. 1999 & 2000) 

The source evades simple classification, although it 
seems to spend most of the time in the very high state. 
Several attempts to categorize the behavior of GRS 
1915+105 have been made in the past. Belloni et al. (2000) 
defined 12 different states depending on light curve varia- 
tion and hardness colors. One of these states, the so-called 
X-state is characterized by a lack of obvious variations in 
the light curve and spectrum and is associated with con- 
tinuous radio emission of differing strength, ^-states corre- 
spond to the low/hard state of GRS 1915+105, exhibiting 
relatively low flux from the accretion disk and showing a 
hard power law tail. Theses states resemble long variants 
of the lulls in /3-states, when part of the inner accretion 
disk is proposed to be absent. Depending on the strength 
of the radio emission, differing phase lag and power spec- 
trum behavior are seen. During the so-called radio quiet 
X-states the phase lag of the hard X-ray photons to the soft 
X-ray photons is partly negative, whereas it is always pos- 
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itive when the radio emission is strong (Muno et al. 2001). 
Also, the frequency of the 0.5-10 Hz QPO decreases with 
increasing radio flux. 

In this paper we present a comprehensive study of the 
X-ray spectral behavior of GRS 1915+105 in the x-state 
with over nearly four years of observation with RXTE. 
The variation of the spectral properties is then compared 
to those of the radio emission. Valuable information about 
the disk geometry and properties of the compact object is 
derived. 

2. Data selection and analysis 

2.1. RXTE 

For our investigation of GRS 1915+105 we used public 
RXTE data from November 1996 to September 2000 pro- 
vided by HE AS ARC. We selected Proportional Counter 
Array (PCA) and High Energy X-Ray Timing Experiment 
(HEXTE) data of 139 x-state observations from 89 differ- 
ent days. The selection of the datasets was based on the 
X-states defined by Belloni et al. (2000) and on PCA light 
curves provided by E. H. MorganQ Data with more than 
30 min away from the South Atlantic Anomaly and with- 
out X1908+075 in the HEXTE background| were selected. 
The lack of obvious variability in the light curve and its 
spectral hardness allows long continuous exposure times 
and therefore high signal-to-noise ratios. Typically, the 
spectra during this state show blackbody emission aris- 
ing from an accretion disk and a dominating power law 
hard energy tail. 

We reduced the RXTE data using the standard re- 
duction script REX included in the HEAsoft5.04 package. 
We restricted the analysis to Standard 2 binned data of all 
layers of PCUO of the PCA from 3-25 keV and HEXTE 
cluster from 20-190 keV only. We used PCARSP 7.10 
to produce a particular response matrix for each PCUO 
dataset. In order to account for part of the uncertainties 
in the PCUO instrument we added a systematic error of 
1 % as recommended (Remillardj^ . 

The X-ray spectral fitting was done using XSPEC 11.0 
(Arnaud 1996). A consistent model should fit all RXTE 
X-state data. Therefore we tested several X-ray radia- 
tion models and finally selected a model consisting of 
(i) photoelectric absorption (WABS; Balucinska-Church & 
McGammon 1992), (ii) a spectrum from an accretion disk 
consisting of multiple blackbody components (DISKBB) 
and (iii) a power law spectrum reflected from an ionized 
relativistic accretion disk (REFSCH; Fabian etal. 1989, 
Magdziarz & Zdziarski 1995). Several attempts in the 
past to fit GRS 1915+105 spectra have shown complicated 
residuals in the soft X-ray band suggesting the existence 
of emission and/or absorption features near 6.4 keV, the 
energy of the Fe Ka line (Kotani etal. 2000). Because of 

^ http:/ /xte.mit.edu/'~ehm/1915_frames.html 
^ http:/ /mamacass.ucsd.edu:8080/cgi- 
bin /HEXTErock.html 
^ http:/ /Iheawww. gsfc.nasa.gov/users/keith/ronr.txt 



the low spectral resolution of '^IkeV of the PCA at this 
energy an additional line fit gives no meaningful results. 
Therefore we ignored all energy bins from 4.5 to 8.5 keV 
while fitting our model. We fixed the hydrogen column 
density at Nh=5-10^^ cm~'^ as determined by Greiner 
et al. (1994) with ROSAT. It has been shown that the disk 
blackbody + power law assumption strongly overestimates 
the flux at lower energies compared to the thermal and 
non-thermal comptonization models used by Vilhu et al. 
(2001) and Zdziarski etal. (2001), respectively. This ex- 
plains the smaller Nh (2-3T0^^ cm~^) values they found 
for x-state observations. However, because we ignored the 
energy bins from 4.5 to 8.5 keV the amount of data bins 
needed to adjust the hydrogen column density was to small 
to let Nh be a free parameter. 

We fixed 9 of the combined 16 model parameters 
(Table leaving free the accretion disk temperature, Thh, 
and relative normalization, Khi, oc , the power law pho- 
ton index, F, with N ocE~'" and relative normalization, 
Kpo (photons/keV/cm^/s at IkeV), the reflection index, 
R, the ionization parameter, ^, and a factor to account 
for the relative normalization between PCUO and HEXTE 
cluster 0. If not stated otherwise, we plot Ict errors for each 
parameter of interest. 

Table 1. Fixed parameters of the spectral model. 



fix. parameter value 

"/Vh 54(Pcm^ 

cutoff energy no cutoff 

redshift 

Z>2 element abundances 1 

iron abundance to abundance above 1 

inclination angle 70° 

power law index for reflection emissivity -2 

inner disk radius 6GM/c'^ 

outer disk radius 1000 GM/c^ 



2.2. GBI and RT 

The appearance of steady radio emission in X"States 
suggested us to search for correlations between RXTE 
data and radio data at 15 GHz {Ryle Telescope = RT) 
and 2.25GHz {Green Bank Interferometer = GBI). GRS 
1915+105 shows variability at all frequencies on time 
scales of seconds to hours. For a useful statement in x~ 
states a suitably selection of corresponding datasets is 
therefore required. 

According to the general interpretation that the radio 
emission is synchrotron emission from ejected plasma in 
sporadic or continuous jets (Fender et al. 1995) the radio 
flux should reach a maximum 15 min after the actual ejec- 
tion (Mirabel et al. 1997) and therefore after a possible 
determining X-ray event. 

Fig. |] shows the 15 GHz radio flux, Fji, from the RT 
for JD 2450898-2450913 together with two RXTE obser- 
vations. The iJT" observed GRS 1915+105 several times a 
day with five minute exposures. 
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Because the radio exposure is much shorter than the 
X-ray exposure ('^hour) it is important to select radio 
data simultaneous with the RXTE observations. 37 of the 
139 analyzed RXTE observations have simultaneous RT 
data and 9 have simultaneous GBI data. 

Still, the selection of simultaneous radio observations 
is non-trivial. Occasionally the radio emission varies also 
during a single x-state X-ray observation, whereas no vari- 
ability is seen in X-ray count rate and hardness ratio (5.2- 
60keV/2-5.2keV) (Fig. |). But the variation of Fr dur- 
ing a x-state observation is negligible compared to the 
uncertainties of the individual radio measurements and 
to the variation between different RXTE observations. 
Therefore, the radio fluxes were averaged for each RXTE 
observation. 
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Fig. 1. Variation of the 15 GHz flux (crosses) between JD 
2450898 and JD 2450913. The dotted Hues mark the time 
of RXTE observations from 29.03.1998 (JD -2450902) 
and 06.04.1998 (JD -2450910). GRS 1915-fl05 showed 
low radio emission until JD 2450900.5. The source was 
dominated by strong, variable radio emission following JD 
2450905. No statement can be made about the time in 
between, because of the lack of radio observations dur- 
ing the RXTE observation from 29.03.1998. Also for the 
RXTE observation from 06.04.1998 the radio flux is un- 
certain because of the strong variation before and after 
the X-ray exposure. 



3. Results of the analysis 

The spectra of all 139 datasets have been fitted with the 
above-described DISKBB-f-REFSCH model. Two example 
model spectra are shown in Fig. ||. These spectra have the 
extreme values of F and i?, revealing the spectral differ- 
ences between different x-state observations. 

The best fit parameters for the x-states of GRS 
1915-1-105 are shown in Fig. and Table H, respectively. 
The upper two panels in Fig. Hshow the strong variability 
of the source in the 1.5-12 keV {RXTE All Sky Monitor 
(ASM)) and 2.25 GHz {GET) bands from JD 2450300- 
2451900. Irregular outburst and relatively quiet phases al- 
ternate in X-ray and radio without any obvious coupling. 

The lower six panels of Fig. ^ show the PCUO count 
rate and the fit parameters of the x-state observations. 
The red and black data points mark two groups of obser- 
vations with different T{Kpo) behavior (see Fig. |[). 



T 1 1 1 1 1 n 1 1 1 1 1 1 1 1 1 1 1 r 




20 40 

Time (min) 



Fig. 2. 15 GHz light curve {RT, upper pannel), 2-60 keV 
PGA light curve (middle pannel) and hardness ratio 
\^-b2keV (lo'^er pannel) of the x-state from 15.09.1998 
(JD 2451071.9). While the radio light curve shows struc- 
tured variability, none is seen in X-ray intensity and spec- 
trum. The large black dot (marked for clarity with the 
dotted lines) in the upper panel represents the averaged 
radio flux over the RXTE observation used for further 
analysis. 



Between the x-states, the PCUO count rate varies, with 
some exceptions, by a factor of 2 at most. Only during the 
long continuous x-state at JD —2450400-2450600 and at 
JD —2451750 recurrent variability in the X-ray count rate 
is observed. The model parameters of DISKBB-I-REFSGH 
are variable. The power law slope, F, varies between 2.4 
and 3.5 with a long-term periodicity of —590 days (Rau 
& Greiner 2002, in preparation). No correlation of F and 
PGUO count rate is seen, except around JD —2450500 and 
JD —2451750. The power law normalization behaves sim- 
ilarly to the slope. 

The reflection component, i?, is variable between dif- 
ferent x-state observations and shows a long-term vari- 
ability similar to the power law slope. It varies between 
and 10 (upper limit of our model) with rather large un- 
certainties when R >4. Except for five observations when 
i? <1, this rules out an isotropically sandwiching corona 
above the entire accretion disk. 

The significance of the accretion disk component varies 
through the x-states. For some observations the disk is 
more or less absent (e. g. JD 50400-50600; relatively low 
disk temperature and large inner disk radius (small nor- 
malization)), for other observations the disk component 
provides a non-negligible contribution to the X-ray flux 
(e. g. JD 51500-51600; high temperature and small inner 
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Fig. 3. Combined PCUO and HEXTE cluster spectra 
of the x-states from 09.10.1997 (down, artificially offset 
for better discrimination) and 30.03.2000 (top). Theses 
spectra show extreme power law slope (r=3.41±0.06 
(09.10.1997), r=2.65±0.05 (30.03.2000) and reflection 
(i?=7.43±1.12 (09.10.1997), i?=0.35±0.23 (30.03.2000) 
behavior. Note the obvious reflection hump for the 
09.10.1997 spectra at ~10keV. 

radius). In order to be consistent and to fit all ^-state ob- 
servations with the same model, we included the DISKBB 
component in all of our fits, although it could be excluded 
for several observations. 



Obs-ID 


K 

[Pho/keV/cmVs] 


^(10) 


c(ii) 


^2(12) 


K-02-01 
K-06-00 
K-07-00 


43.51^;^ 

le.slji 


4999 3521 

50001:^421 


0.82 
0.78 
0.79 


0.69 
1.00 
1.01 


P-28-00 
P-28-01 
P-28-02 


04 1+0.7 

29.8tli 
27.2tl:l 


2902t^2^«o« 

q/i+2900 


0.82 
0.81 
0.52 


0.94 
0.84 
0.77 



The accretion disk component shows a variable disk 
temperature of 1-4 keV. Sometimes, large uncertainties 
due to the small contribution of the DISKBB component 
to the total flux are seen. The inner disk radius, which 
can be determined from the disk normalization, varies 
between 1 and 20 km. For a non-rotating black hole of 
mass 14 Mq (as measured for GRS 1915-1-105, Greiner 
etal. 2001b) the Schwarzschild radius is ~40km. It is 
known that the DISKBB model underestimates the in- 
ner disk radius by a factor of 1.7-3 due to Doppler blur- 
ring and gravitational redshift (Merloni etal. 2000). Also 
the neglect of comptonization in the surface layers of the 
disk leads to unphysical values when using then DISKBB 
model (Zdziarski etal. 2001). But even a maximally ro- 
tating black hole (where the inner disk radius reaches the 
Schwarzschild radius) cannot account for the majority of 



the parameter values. This problem has to be kept in mind 
when discussing the absolute values of the parameters. 

Another free physical parameter of the REFSCH 
model is the ionization parameter, ^. It has huge uncer- 
tainties because no Fe Ka line could be fitted but has a 
negligible infiuence on the hard X-ray continuum and our 
model parameters. Therefore we will not plot or discuss 
^ further. Theoretically, the reflection in the hard spec- 
trum should depend on the ionization parameter because 
higher ionization means lower absorption and therefore 
higher reflection probability. But the REFSCH model in- 
cludes a simple 1-ionization zone model only, which is very 
unlikely to be present in the disk and does not show any 
dependence of reflection and ionization at all. 

3.1. The power law component 

The model fits reveal an increasing power law normaliza- 
tion, Kpo, with a steepening power law component (Fig.|5|) 
suggesting a pivoting behavior. Two branches with differ- 
ent slopes are seen in the correlation. No correlation of the 
power law slope with the X-ray count rate in ASM and/or 
PCU is observed. 




Fig. 5. Power law slope, F, as a function of power law nor- 
malization, Kpo. Each RXTE observation is represented 
by one data point. The upper branch presents the red 
points in Fig. ^ and the lower branch the black points, re- 
spectively. The dotted lines represent the fltted correlation 
functions (see text). 

The strength of the correlation can be tested using a 
Spearman rank-order correlation test (Press etal. 1992). 
Both branches show strong correlations (steeper: r5=0.85, 
flatter: rs=0.87). (Note, this statistical method does not 
take into account the particular uncertainties of the data 
points.) 

The best descriptions of the correlations are functions 
of the type 

T = u- logK + V (1) 
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Table 2. Fit results of the RXTE spectra using the DISKBB+REFSCH model (The complete table is available 
at CDS). (1): ID of observation (1=10402-01, J=20187-02, K=20402-01, L=30182-01, M=30402-01, N=30703-01, 
0=40703-01, P=50703-01), (2): exposure time of observation, (3): 3-20 keV PCUO count rate, (4): 20-190 keV HEXTE 
cluster count rate, (5): accretion disk temperature, (6): accretion disk normalization {Kbb = ( d/io* fcpc )^ cos9), (7): 
power law slope, (8): reflection amplitude, (9): power law normalization (photons/keV/cm^/s at 1 keV), (10): ionization 
parameter, (11): factor for normalization of PCUO and HEXTE cluster 0, (12): reduced . All errors are la for each 
parameter of interest. 

Obs-ID^^) GD JD Exposure^^' PCU(F^ HEXTECF^ R**' 
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K-02-01 
K-06-00 
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14.11.1996 
11.12.1996 
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with M=0.57±0.02 and ?;=2.0±0.1 for the lower branch 
and M=0.81±0.07 and t;=2.0±0.1 for the upper branch. 

The pivoting has been tested in more detail by 
schematically plotting the obtained power law spectra. 
The two different branches consequently show different 
pivoting behavior (Fig. H) . Whereas the upper branch from 
Fig. ^ shows a pivoting energy of 4-8 keV (gray) , the lower 
branch pivots at around 20-30 keV (black). 




Energy (keV) 

Fig. 6. Schematic power law spectra for the upper (top) 
and lower (bottom) branch of Fig.|[ the red and black 
points from Fig. ^, respectively, artificially offset by a fac- 
tor of 100. The pivoting happens at different energies (4- 
8keV & 20-30 keV). 



3.2. Reflection 

Due to the known uncertainties of the PCA below 6keV 
(Jahoda et al. 1996) and the low energy resolution, strin- 
gent conclusions about the behavior of an iron line in the 
RXTE spectra at 6.4-7 keV not can be drawn. For the 
analysis of the reflection, therefore, only the continuum 



radiation can be used, manifesting itself in the reflection 
hump from 10-30 keV. 




1.5 2 2.5 3 3.5 

r 



Fig. 7. i?(r)-correlation for GRS 1915+105 (crosses) and 
three X-ray binaries (Cyg X-1 = circles, GX 339-4 = tri- 
angles and GS 1354-644 = squares; from Gilfanov et al. 
2000). The solid line represents the best fltting model 
function and the dotted lines the upper and lower lim- 
its. Overplotted are the 1, 2 and 3(t confidence contours 
for three different observations (from left: N-23-01, K-49- 
01, O-40-03). Note the data points (JD 2451497-2451513) 
at r '^3.3-3.5 and R ^0-2 behaving remarkably different. 
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The reflection amplitude in GRS 1915+105 as a func- 
tion of r is shown in Fig. ^ together with that of three 
black hole candidates (Cyg X-1, GX 339-4 and GS 1354- 
644; Gilfanov etal. 2000). Although R shows large uncer- 
tainties for the x-states of GRS 1915+105, a similar cor- 
relation is seen as in the other X-ray binaries. The steeper 
the power law component, the higher is the reflected frac- 
tion of photons. Fig. |^ contains the confidence intervalls 
of R for three observations to clarify the influence of the 
model on the i?(r)-correlation. Although the contours are 
elongated similiar to the correlation (higher F has higher 
R), they are intrinsically steeper compared to the overall 
correlation. Thus, the existence of the i?(F)-correlation in 
the x-states of GRS 1915+105 is no artifact of the model. 

A group of observations at high F (3.3-3.5) and small 
R (0-2) behaves remarkably differently. These points be- 
long to five datasets from JD 2451497 and JD 2451513 
showing very high PCUO count rates (1300-2100 cts/s) 
and are short duration x-states between different high 
variability states. 

The correlation of R and F is tested using a Spearman 
rank-order correlation test. The correlation is distinct 
(rs=0.61) but not strong. 

For a quantitative description of the correlation, a phe- 
nomenological function was fitted to the data. The best- 
fitting model is a power law 

i? = u-F", (2) 

with M=0.05±0.007 and w=3.7±0.4 (solid line in 
Fig.0). 

3.3. X-ray-Radio-Correlation 

No obvious correlation between the soft X-ray component 
and the radio emission is seen. Neither the soft X-ray flux, 
the temperature of the accretion disk nor the inner disk 
radius show a correlation with the 15 GHz radio fiux from 
i?T and/or GBI . 

An unexpected result is found when plotting F vs. Fjf 
(Fig. H). The power law slope correlates positively with 
the radio fiux at 2.25 GHz and 15 GHz. Observations with 
high radio emission show a softer X-ray spectrum (steeper 
power law component). Note that no bimodality in the 
radio emission exists. Instead, a continuous spread is ob- 
served. A strict separation of radio loud and radio quiet x- 
states, as done before (e. g. Muno et al. 2001, Trudolyubov 
2001) seems therefore unsubstantiated. 

The correlation of F and F/j is most obvious for si- 
multaneous (Ai=Ohrs) observations. Datasets with radio 
observations ±5 hrs offset also show the correlation (Fig. || 
upper panel) but with significant scattering. Usually, no 
statement about the variability state before and after the 
RXTE exposure can be made, because the source may 
have had several state alterations. Thus, non-simultaneity 
is the likely reason for the strong scattering of the corre- 
lation for Ai=5hrs in comparison to Ai=Ohrs. 

To determine the strength of the correlation a 
Spearman rank-order correlation test was made for the 
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Fig. 8. F(FR,)-correlation for GBI (left) and RT (right). 
Upper panels: correlation for Ai=5 hrs (time between 
RXTE and radio observation). As for GBI and for RT 
the correlation is visible but the data points show a large 
spread. Lower panel: correlation for Ai^O hrs. The dotted 
lines represent the correlation functions. 

RT data with Af;=Oh. The correlation is strong (r^ = 
0.83). No test was made for GBI due to the small number 
of simultaneous GBI/ RXTE datasets. 

The next step was to fit a phenomenological function 
(dotted line in Fig. ||) to the correlation. The best model 
functions are of the form 



T = u ■ logFij 



V, 



(3) 

RT and 



with u=0.54±0.02 and w=3.94±0.03 for 
M=0.32±0.05 and u=3.3±0.1 using GB/data. 

Besides the FR(F)-correlation we also investigated the 
relation between Fj^ and the X-ray flux of the power law 
component. This introduces some freedom as to at which 
lower energy the power law is "chosen" . Fig. ^ shows the 
X-ray flux from 20-200 keV and 1-200 keV, respectively, 
over the radio flux. 



20-200 keV 



1-200 keV 



510-i- 
410-i- 

E 

\. 310-i- 
w 

X 
[x. 

110-S- 
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0.001 0.01 0.1 0.001 0.01 0.1 

15 GHz Flux (Jy) 

Fig. 9. 20-200 keV (left) and 1-200 keV (right) X-ray flux 
shown over the 15 GHz radio flux (RT). The 20-200 keV 
flux has a negative correlation with the radio flux, whereas 
the 1-200 keV flux does not show such a correlation. 

With increasing radio emission, the X-ray flux in the 
power law component between 20-200 keV decreases (ne- 



A. Rau & J. Greiner: Comptonization and Reflection in GRS 1915+105 



7 



glecting reflection). The total 1-200 keV X-ray flux in the 
power law component is more or less constant. No correla- 
tion with the radio flux is seen. This is because the flux in 
the 1-20 keV range by far dominates, thus washing out the 
correlation. Muno etal. (1999) did not find a correlation 
of Fx (50-100 keV) with the 15 GHz radio flux. But they 
used a model consisting only of DISKBB and a broken 
power law (BKNPO) for the X-ray spectra and RT data 
with±12hrs offset. 

A Spearman rank-order correlation test for Fx(20- 
200) gives r's=-0.75, thus showing a strong anti- 
correlation. For (1-200 keV) the test results in rs=-0.09 
and ratifies the lack of a correlation between the 1-200 keV 
non-thermal X-ray flux and F/j. 



tail up to IMeV without a cutoff (lyudin 2000, Zdziarski 
et al. 2001). But as several contemporaneous RXTE obser- 
vations show (which cover ^3% of the OSSE exposure), 
during the observation GRS 1915-f 105 went through sev- 
eral x/a state transitions and the PGA photon count rate 
from 2-40 keV varied between 5 and 25 kcts/s. Thus, the 
OSSE spectrum is the sum of the spectra of the x ^-nd 
a (and possible other) states, with an unkown fraction 
coming from the x-state. It is therefore premature to un- 
equivocally relate the hard MeV tail to the x-state. 

4.1.1. Bulk Motion Comptonization 



Table 3. Observed correlation of the x-states of GRS 
1915+105. 



correlation between: 



type 



ITT 



ref. 



r Kj,o 

r Fr 

V R 

Fx (20-200 keV) Fr 

J'QPO (0.5-10 Hz) Fr 

i^QPO (0.5-10 Hz) Tbb 



+ 
+ 
+ 



+ 



Muno et al. 2001 
Muno et al. 1999 



(1) : direction of correlation ("+" = positive, "-" = negative), 

(2) : frequency of the 0.5-10 Hz QPO. 



4. Discussion 

Table ^ summarizes the correlations found here plus 2 rel- 
evant correlations found by Muno etal. (1999, 2000). In 
the following we will discuss implications of these correla- 
tions. 

4.1. Comptonization 

The origin of the hard power law component in AGN 
and X-ray binaries is still under discussion. Three main 
models exist, where the hard photons originate by inverse 
Gompton scattering of soft disk photons on hot electrons. 
The main difference between these models is the distribu- 
tion of the electrons. They can be thermal (Maxwellian), 
non-thermal (power law like) or free falling from the last 
stable orbit onto the event horizon of the black hole. The 
high energy spectrum presents an important test for the 
distinction between these models. Therefore, the behavior 
of the power law component in GRS 1915+105 is crucial 
for the understanding of the electron distribution and for 
the geometry of the system. 

It is still unclear whether or not the hard spectrum of 
GRS 1915+105 in the x-state extends with or without cut- 
off up to MeV energies. Due to the rapid variability of the 
source and the long required exposure times it is nearly 
impossible to get sufficiently accurate 7-spectra with the 
present high energy satellites. For example, the OSSE 
spectrum from 14th-20th of May 1997 reveals a power law 




Fig. 10. Bulk motion geometry. Electrons become spheri- 
cal free-falling between the last stable orbit and the event 
horizon due to shocking. The photons arc comptonizcd on 
these accelerated electrons. 

One model for the comptonization is the bulk motion 
comptonization (Blandford & Payne 1981, Ghakrabarti & 
Titarchuk 1995) (Fig. The accretion stream passes 
the transition radius, r^s, at the last stable orbit around 
the black hole and gets shocked. The matter falls spherical 
onto the event horizon. In front of the shock, the stream 
is similar to an optically thick, geometrically thin accre- 
tion disk (as required by the strong observed blackbody 
component). 

The free falling electrons are accelerated up to the 
speed of light and inverse Gompton scattering of soft pho- 
tons on the electrons provides the observed power law in 
the hard spectrum. The model predicts a cutoff depend- 
ing on the mass accretion rate (Ebisawa et al. 1996) due 
to inverse Gompton scattering and Gompton recoil. As 
stated above, the existence/absence of a cutoff during the 
X-states is still not secured, wherefore not statement about 
the model can be made from that issue. 

Two major problems exist for the description of x~ 
state data by the BMG model. The RXTE spectra of GRS 
1915+105 partly show a strong reflection component. This 
is not compatible with the geometry of a thin disk out- 
side a central, spherical accretion stream. The accretion 
stream is accelerated away from the inner disk edge onto 
the event horizon of the black hole. Therefore only a small 
part of the soft photons are comptonized back into the disk 
plane at r > rts- In conflict to the observations, no strong 
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reflection component can emerge. A possible explanation 
for the observed reflection can then be a partly covering of 
the disk by an absorbing medium with tt '^3 (Zdziarski 
2000). Recent Chandra X-ray spectra indeed show a high 
X-ray column density and abundance excesses for Si and 
Fe, which may be related to material that is associated 
with the immediate environment of the source (Lee et al. 
2002). 

The other problem for the BMC model is the time lag 
of hard and soft X-ray photons. Muno etal. (2001) found 
a phase lag in the 5 Hz QPO (which are connected to the 
hard X-ray photons) of 0.5. This corresponds to a delay of 
0.1 s of the hard photons to the soft photons. The expected 
delay for scattering on a convergent electron stream inside 
rts is ^0.2 ms, much smaller than observed. The phase lag 
is not constant in time and frequency and seems to vary 
with the radio flux. 

It is obvious that the BMC model not can explain 
the hard power law component in the X"States of GRS 
1915+105. There are indications that the BMC model 
fits other variability states. Shrader and Titarchuk (1998) 
found good spectral agreement of the BMC model and 
RXTE observations in the p-state. Note, however, that 
they completely neglected the time delay of hard and soft 
photons and the strong variability of the p-statc for their 
X-ray spectral modelling. 

4.1.2. Thermal Corona 

Another possible model for the origin of the hard X-ray 
component is the disk corona geometry (Fig. |ll|) (e.g. 
Haardt & Maraschi 1993, Svensson & Zdziarski 1994). 
Soft disk photons are inverse Compton scattered on hot 
(>20keV) electrons located above the disk. 



likely due to thermal comptonization. The predicted cutoff 
in the hard X-ray spectrum is seen in Cyg X-1 (Gierliiiski 
et al. 1997), in the LMXB GX 339-4 (Zdziarski et al. 1998) 
and some other X-ray transients (Grove etal. 1998). 

Also the high/soft state in X-ray binaries can be ex- 
plained by thermal comptonization but requires high elec- 
tron temperatures in some systems to explain the observed 
unbroken power law up to MeV. In addition tt ^1 is re- 
quired in order to keep the spectrum soft (Zdziarski 2000). 

The power law slope is a measure of the Compton am- 
plification, A{Lsoft + Ldiss)/Lsoft- Ldiss is the dissipated 
energy in the corona and Lsoft the energy of the incoming 
photons. Table ^ shows how A depends on F with 

T^l{A-l)-' (4) 

(Beloborodov 1999) and 5—1/Q for galactic black holes. 
The ratio of the dissipated energy in the corona and the 
energy of the intercepted soft photons gets smaller with a 
steeper spectrum. For the highest F only ~9% of the soft 
photon energy is dissipated into the corona. 

Table 4. Fraction of the dissipated energy in the corona, 
Ldiss, on the energy of the incoming photons, Compton 
parameter, y = AT^kTe/ (nieC^) {tt is the Thomson 
depth), and Compton amplification. A, for three different 
power law slopes. 



F 


Ldiss 


y 


A 


2.4 


84% 


1.33 


1.84 


3.0 


22% 


3.65 


1.22 


3.5 


9% 


7.30 


1.09 




Fig. 11. Disk corona geometry. The disk photons are 
comptonized in an electron distribution located above the 
accretion disk. (Note, the shape and size of the electron 
distribution in this plot has no real physical or observa- 
tional meaning.) Part of the photons are scattered back 
into the disk plane and are reflected or reprocessed. The 
photons reaching the detector are compound of the soft 
disk photons, the comptonized and the reflected photons. 



Simultaneous CC?i?0/OSSE and X-ray observations of 
X-ray binaries with black holes showed that the high en- 
ergy continuum in the low/hard state originates most 



The minimum energy which a photon receives when 
passing a thermal electron distribution depends on the 
electron temperature fcTg 

— = 2 • W 

The result is a lack of low energy photons when kTg. be- 
comes large. For fcTe>200keV the discrete orders (mainly 
the first) of the Compton scattering become visible as 
an additional peak overlaying the soft disk component. 
The high disk temperatures (3-4 keV) in GRS 1915+105 
provide a photon deficit below 6-8 keV. This makes obvi- 
ous, as already mentioned above, that a simple power law 
model is insufficent for the description of the comptonized 
photons. The lack of photons at low energies affects F, 
Kpo and the reflection component in the model. 

4.1.3. Hybrid Corona 

The radiation processes and geometry can also be de- 
scribed by comptonization of soft disk photons on a hybrid 
(thermal and non-thermal) electron distribution above an 
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optical thick accretion disk (as in Fig. |Tl|)(Coppi 1992). 
Part of the electrons in the thermal corona are acceler- 
ated, possibly in reconnection events. The disk photons 
are comptonized on these electrons forming the observed 
power law component in the hard X-ray spectrum. The 
energy of the non-thermal electrons is partially trans- 
fered to the thermal electrons due to Coulomb scattering. 
This leads to heating of the thermal electrons above the 
Compton temperature. Then the thermal electrons play 
an important part in the up-scattering of the soft disk 
photons, too. 

The compactness (ratio of luminosity and size) is an 
important parameter of the coronal plasma (Coppi 1999). 
A large compactness leads to electron positron pair cre- 
ation due to photon photon collisions resulting in a pair 
annihilation line at 511 keV in the spectra. When the 
compactness is low, the loss of energy of the electrons 
due to Coulomb scattering dominates the loss due to in- 
verse Compton scattering. The ratio of thermal to non- 
thermal electrons increases and the plasma becomes ther- 
mally dominated. Together with the resulting break in the 
distribution of the non-thermal electrons this results in a 
cutoff of the photon spectrum. 

The rate at which non-thermal electrons appear in the 
hybrid electron distribution can be written as a power law 
(Gierhhski ct al. 1999) 
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The investigation of the RXTE data of GRS 1915-^105 
with the DISKBB+REFSCH model docs not allow a defi- 
nite conclusion about the origin of the hard spectral com- 
ponent. Only the BMC model is ruled out for the x-states. 
Though a clear distinction between a thermal and a hy- 
brid not can be proposed, it seems clear that the thermal 
electrons dominate. 

It has been shown above (Fig. P) that two different 
states of comptonization with two different pivoting en- 
ergies are observed. Very recently (Zdziarski ct al. 2002) 
pivoting was found also in Cyg X-1. A varying amount 
of soft seed photons undergoing the comptonization pro- 
cess probably accounts for the pivoting of the spectrum. 
One possible explanation is a constant disk blackbody and 
an expanding/contracting hot coronal plasma. A change 
in the size of the corona leads to a change in the frac- 
tion of intercepted soft photons. Assuming constant op- 
tical plasma depth, tt, if the soft luminosity is high, the 
electrons in the corona are cooled efficiently and the spec- 
trum of the comptonized photons is soft. On the other 
hand, if the soft luminosity decreases, the temperature 
of the corona increases and the spectrum hardens). The 
ocurrance of two differing pivoting branches suggests the 
existence of two states with different tt and/or different 
electron composition (thermal/non-thermal). 



Here, 7 is the Lorentz factor of the non-thermal electrons 
and Tin the power law slope of the injected soft photon 
distribution, which influences the slope of the power law 
spectrum as 

^ 2(F - 1). (7) 

The observed power law slope F ~2.5 (F '^3.5) implies 
Nnt 7^"^ (oc 7^^). Therefore the fraction of non-thermal 
electrons decreases with steepening power law. 

In GRS 1915-1-105, the observed continuous distribu- 
tion of F between 2.4 and 3.5 is evidence for a variability in 
the power law component and thus in the composition and 
locus of the electron distribution. From the observation of 
the reflection and the 0.5-10 Hz QPO it is conspicuous 
that the electron distribution must have a small (relative 
to the system) spatial size. As a small compactness, e. g. 
large plasma volume, results in a spectral cutoff, for the 
X-states of GRS 1915-1-105 a hybrid electron distribution 
not can be ruled out. It is difhcult to distinguish spectro- 
scopically between a hybrid and a thermal plasma if the 
hybrid plasma is thermally dominated. The only possibil- 
ity is the search for the predicted photon excess at high 
energies and for the annihilation line (Coppi 1999). For 
a better understanding, enhanced high energy detectors, 
such as on INTEGRAL, are required. The analysis of such 
spectra also needs further developed physical models, in- 
stead of a simple power law. 



4.2. Reflection 

The reflection of X-rays on an accretion disk manifests it- 
self in two specific spectral features: (i) the characteristic 
emission line spectrum (mainly the Ka lines of the most 
common metals) with the 6.4 keV iron emission lines as 
the strongest and (ii) a characteristic hump arises above 
lOkeV because of the energy dependence of the cross sec- 
tion of absorption and Compton scattering. 

Observations with ASCA in different variability states 
indicate a variable iron absorption or emission in GRS 
1915-1-105 (Ebisawa etal. 1997). The iron emission at 
6.4 keV dominates the spectra from 25.10.1996 and 
25.04.1997 whereas a distinct absorption at 7keV and no 
emission was seen on 27.09.1994 and 26.04.1995. During 
the 25.04.1997 observation GRS 1915-^105 switched be- 
tween X- and a-states (based on RXTE data). Therefore, 
it is not clear whether the emission line originates during 
the X- or during the a-interval. No statement can be made 
about the states during the other ASCA observations due 
to lack of corresponding RXTE observations. 

Recent Chandra data of GRS 1915-1-105 obtained in 
the low hard state revealed neutral K absorption edges, 
ionized resonance absorption from Fe (XXV, XXVI) and 
possible emission from neutral Fe Ka and ionized Fe XXV 
(Lee et al. 2002), suggesting conditions favorable for reflec- 
tion. 
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4.2.1. Reflection in LMXB and AGN 

Investigations of Seyfert galaxies and X-ray binaries in the 
low/hard state show a correlation between the slope of the 
hard X-ray component and the reflection. Incoming radi- 
ation with a steeper slope is more reflected than radiation 
with a flatter slope. 

The reflection is stronger in black hole candidates than 
in Seyfert galaxies for a given T (Zdziarski 1999). This can 
be explained using an optical thick accretion disk. The dif- 
ferent masses of the black holes imply different maximum 
energies of the blackbody photons (Svensson 1996). The 
hard X-ray photons from the corona can ionize the up- 
per layers of the accretion disk and an ionized layer above 
neutral matter evolves. If F is large the heating of the disk 
is small and the cold layer of the disk lies near the disk 
surface. With decreasing power law slope the temperature 
of the upper layer and therefore the ionization increases, 
the absorption is smaller, more photons are scattered and 
the characteristic reflection hump is suppressed, sufficient 
to explain the observed _R(r)-correlation. 

The correlation in Seyfert galaxies is generally steeper 
than that observed in X-ray binaries (Zdziarski 1999). 
In Seyferts the correlation has been explained with a 
model consisting of a static, thermal corona above a 
neutral reflector (Svensson 1996). The existence of the 
R(r)-correlation implies a feedback where the existence 
of the reflecting matter influences the hardness of the X- 
ray spectrum (Bottcher etal. 1998, Zdziarski etal. 1999). 
Assuming that the cold matter (the accretion disk) emits 
soft photons which become seed photons of the comp- 
tonization, than with increasing solid angle of the reflector 
(here the accretion disk) the flux of soft photons increases 
and the cooling rate of the hot corona increases. For a 
thermal plasma the resulting power law component steep- 
ens with increasing cooling rate. 

On the other hand, models which are based on non- 
thermal electrons should not show a dependence of spec- 
tral hardness and cooling rate (Lightman & Zdziarski 
1987, Zdziarski et al. 1999). Models in which the seed pho- 
tons are intrinsically produced in the hot plasma (e. g. 
synchrotron radiation) not can reproduce the observed 
correlation. The power law slope is then independent of 
the reflection. 

4.2.2. Reflection in GRS 1915-^105 

The x-states in GRS 1915-F105 show a similiar correlation 
as seen in Seyfert galaxies and X-ray binaries (Zdziarski 
etal. 1999). The reflection amplitude is strongly effected 
by the disk model. The higher the disk temperature, the 
larger the reffection amplitude. The disk temperatures in 
X-ray binaries are higher than in Seyfert galaxies. In GRS 
1915-1-105 the temperature of the disk is still higher, ex- 
plaining why a majority of the GRS 1915-1-105 data lies 
above the correlation function (showing slightly higher R 
for a given T) in Fig. 0. 



The reflection amplitude of i? > 1 strongly suggests 
an anisotrophic inverse Compton process. The dominant 
fraction of the up-scattered soft disk photons is directed 
back into the plane of the accretion disk. This produces 
the observed large reflection amplitude in contrast to an 
isotrophic scattering where R < 1. 

Observations of GRS 1915-1-105 show a strongly vary- 
ing r, therefore the corona can hardly be static above 
the reflector. In GRS 1915-1-105, F increases with increas- 
ing radio emission. Higher radio emission may imply more 
outflow away from the disk with f3 = {v/c) >0. Thus, the 
higher the mass outflow, the lower the electron temper- 
ature which is required to produce the observed amount 
of comptonization. On the other hand, with a dynamical, 
thermal corona the spectrum hardens with increasing /3 
due to relativistic aberration, giving a flatter power law 
with higher outflow velocity (Malzac etal. 2001). 

The observation in GRS 1915-1-105 predicts an at least 
partly thermal source of the hard X-ray photons and an 
important contribution of the cooling due to the soft pho- 
tons. This allows both a thermal and a hybrid electron 
distribution. 

It is more difficult to interpret the behavior of observa- 
tions with high F and small R in Fig. ^ Unlike in a thermal 
plasma, the Compton scattering in a non-thermal plasma 
does not depend on the coohng rate of the soft photons, 
instead it depends on the slope of the electron distribution 
(Poutanen & Coppi 1998). A change in the properties of 
the plasma can provide the differing behavior. 

For a better determination of the correlation and for 
conclusions about the distinct system components, more 
complex ionization models are needed. They partly exist 
(Done & Nayakshin 2000) but require too much computing 
power to be useful for the analysis of real spectra at this 
stage. 

4.3. X-ray-Radio-Correlation 

The lack of a correlation of the radio flux with any of 
the disk parameters is somewhat unexpected since the 
(short-duration) low/hard states have earlier been related 
to jet ejections with radio emission by synchrotron radia- 
tion of the ejected plasma (Pooley & Fender 1997). Since 
IR (Eikenberry et al. 2000) and radio observations (Pooley 
& Fender 1997, Fender & Pooley 2000) have led to the con- 
clusion that GRS 1915-1-105 shows jets on various scales, 
or even has a continuous distribution of jet strength, one 
would have expected that also in ^-states the radio flux 
correlates with the disk temperature/emissivity. 

The matter in the sporadic, relativistic jets may orig- 
inate from disruption of the inner part of the accretion 
disk during /3-states (Mirabel etal. 1998). In x-states no 
correlation of the radio emission and the inner disk ra- 
dius is observed. Instead, the hard spectral component is 
correlated with the radio flux. 

Muno etal. (1999) found a positive correlation of the 
accretion disk temperature and the frequency of the 0.5- 
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10 Hz quasi periodic oscillations, and of the radio flux with 
the QPO frequency (Muno etal. 2001). With increasing 
Tbb and decreasing radio emission the QPO frequency in- 
creases. This predicts a negative correlation of Tbb and 
radio flux. The fact that nothing like this is found in 
the present analysis can have several reasons. As men- 
tioned before, Muno et al. (1999) used the standard model 
DISKBB+BKNPO for the analysis of the RXTE spectra. 
This seems not suitable for the x-states of GRS 1915+105. 
Accordingly, the disk temperature and therefore the cor- 
relation of the QPO frequency and the disk temperature 
have to be interpreted carefully. 

It is generally assumed that both the radio jets (Fendt 
& Greiner 2001, Fender 2001) and the hard spectral com- 
ponent originate near the black hole. With increasing out- 
flow of matter (therefore increasing radio emission) the 
interaction with the electron distribution becomes signifi- 
cant. Similar to the interpretation of the pivoting behav- 
ior of the X-ray spectra, an increasing outflow of matter 
implies an increasing size of the scattering medium and 
therefore an increasing amount of intercepted soft seed 
photons. This leads to a lower plasma temperature due 
to cooling and to a softer X-ray spectrum, which is seen 
in the r(_Ffj')-corrclation. Simultaneously, the outflowing 
matter intermingles with the coronal matter and pushes 
it away from the accretion disk. The formerly thermally 
dominated electron distribution may become non-thermal 
dominated. This should result in a shift of the cutoff of 
the power law component in the x-states with higher ra- 
dio emission above the HEXTE range to higher energies. 

When the radio emission is low, the observed power law 
is supposed to originate in the corona above the accretion 
disk. Assuming that the radio quiet state lacks outflow of 
matter or the effect of the outflow on the corona is neg- 
ligible, the hard spectral component possibly originates 
due to comptonization of disk photons in the corona. A 
thermal electron distribution can then explain the X-ray 
spectra of the radio quiet x-states. 

The hard X-ray component in x-states therefore comes 
most likely from soft disk photons which are inversely scat- 
tered on a thermal dominated electron distribution (when 
radio flux is low) or on the base if a continuous jet (when 
radio flux increases). The suggestion of the base of the jet 
as source of the hard X-ray photons was already made by 
Fender (2001). Only tlu^ increasing reflection with increas- 
ing radio flux is still unresolved. 

In conflict with the anti-correlation of Fx (20-200) and 
Fn in GRS 1915-f 105 a linear correlation is observed on 
the LMXB GX 339-4 (Corbel et al. 2000). With increasing 
radio flux the X-ray flux increases. This is true for the hard 
(20-100 keV) and for the soft X-ray flux (2-12 keV) where 
GRS 1915-1-105 shows no correlation. The soft disk compo- 
nent in GX 339-4 is negligible above 2 keV and the power 
law component alone quantifies the flux. Although GRS 
1915-1-105 shows an equal i?(r)-correlation as GX 339-4 
the T{Fi{)- and Fx-(Fi{)-correlation is opposite. That is 
remarkable because the i?(r)-correlation suggests a simi- 
lar structure near the black hole (especially the corona) in 



GRS 1915-hl05 and GX 339-4. No such correlation of Fx 
and Fji is observed in Cyg X-1 (Brocksopp etal. 1999). 

In principle the hard X-ray photons can originate in 
the continuous jet itself. But following the equations of 
Marscher (1983) the hard X-ray photons from self comp- 
tonization of a compact syncrotron jet are negligible with 
regard to the coronal contribution. An analogous result 
was found for Cyg X-1 and GX 339-4. The contribution of 
thermal X-ray bremsstrahlung in the jet is also too small 
to explain the observed hard X-ray luminosities (Memola 
etal. 2002). 



5. Conclusion 

The analysis of 139 EX TE^ observations of GRS 1915-M05 
in the x-state revealed variable components and parame- 
ters of the X-ray spectra. The structure of the hard X-ray 
power law depends e. g. on the radio flux. Further, a two- 
branch correlation of the power law slope and the power 
law normalization was found. Theses branches show dif- 
ferent pivoting behavior. 

The most probable geometry is that of a hot corona 
above the accretion disk. A bulk motion comptonization 
seems to be ruled out for the x-states because of the time 
lags of hard and soft X-ray photons. 

The continuous outflow of matter is not correlated 
with the accretion disk parameters as measured at X-rays. 
Neither the disk temperature nor the inner disk radius are 
connected with the radio flux. Therefore, the outflowing 
matter should not be provided by a disruption of the in- 
ner part of the accretion disk, as thought to be the case 
for the sporadic jets in GRS 1915-1-105, and a positive 
correlation between radio emission and power law slope is 
found. Because both the base of the jet and the corona are 
probably located near the black hole the correlation indi- 
cates an interaction of both structures. Whereas during 
low radio emission the hard X-ray component originates 
due to comptonization on a thermally-dominated corona, 
in radio loud states the comptonization should appear in 
the outflowing matter. 

GRS 1915-1-105 shows a positive i?(r)-correlation as 
seen in other X-ray binaries and AGN. The large reflection 
amplitude suggests a highly anisotrophic inverse Compton 
scattering with the dominant part of the soft photons be- 
ing scattered back into the disk plane. 
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Fig. 4. 1.5-12 keV ASM count rate (1st panel from top) and 2.25 GHz GBI flux (2nd) from November 1996 to 
September 2000. {GBI was off-line some time before JD ^--^2450400 and following JD ~2451660. The lower panels 
show the 3-20 keV PCUO count rate (3rd), the power law slope, F, (4th) and power law normalization, Kpo, (5th), the 



